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In this work we study the propagation of massive Dirac neutrinos in matter with flavor mixing, 
using statistical techniques based on Relativistic Wigner Functions. First, we consider neutrinos in 
equilibrium within the Hartree approximation, and obtain the corresponding relativistic dispersion 
relations and effective masses. After this, we analyze the same system out of equilibrium. We verify 
that, under the appropiate physical conditions, the well known equations for the MSW effect are 
recovered. The techniques we used here appear as an alternative to describe neutrino properties 
and transport equations in a consistent way. 

I. INTRODUCTION 

Neutrino propagation in dense media becomes an important issue in some astrophysical scenarios, such as super- 
novae, neutron stars during the Kelvin-Helmholtz epoch and the solar neutrino problem 0, [Q]. The first two cases 
correspond to compact stars, where densities a few times the nuclear saturation density are reached. To describe neu- 
trino propagation in the dense core of such compact objects, aside of production and absorption ratesjj], one is mostly 
interested in the neutrino cross section with the surrounding matter, which in turn is related to the imaginary part 
of the forward amplitude, due to the optical theorem. In the case of solar neutrinos (and also for the 'atmosphere' 
of compact stars), the possibility of neutrino oscillations and flavor conversion appears. In the standard picture, 
these oscillations are associated to the real part of the forward scattering amplitude of massive neutrinos with the 
background. 

In this paper we will be concerned about the latter topic, i. e. matter effects on neutrino masses and neutrino 
oscillations. Since the seminal papers by Wolfenstein || and by Mikheyev and Smirnov ||, 0, there have 
been many papers which have explored the physics of in-medium neutrinos, specially in connection with the solar 
neutrino problem (for a review, see for example ||] and references Q, Q). Different techniques have been used to 
approach this problem. The simplest one consists in describing the matter effect by an effective potential (as in the 
Wolfcnstein's papers) which will be added to the mass matrix to give an effective Hamiltonian in the Schrodinger's 
equation. Although this approach will suffice for most purposes, it is clear that it does not describe matter effects 
in a covariant way. In order to obtain covariant equations, one has to obtain the corresponding dispersion relations, 
which are the in-medium analogous to the simple mass-shell condition p 2 = m 2 of free particles. Dispersion relations 
of neutrinos in different backgrounds, taking into account the mixing among generations, have been investigated by 
Notzold and Raffelt jl0|]. In that work, dispersion relations appear as the poles of the neutrino propagator (Green 
function), evaluated at the one- loop approximation. 

On the other hand, one has astrophysical situations in which contribution of neutrinos to macroscopic magnitudes, 
such as the energy, pressure, etc, ... becomes important (this is the case in a supernova collapse ||, or in the Early 
Universe). In this case, one has to introduce a distribution function for neutrinos to describe the number of neutrinos 
having a given momentum. It is then desirable to develop evolution equations for these functions in the case where 
neutrino oscillations are present, which implies that distribution functions become non-diagonal matrices in flavor 



space. Several works have implemented this in different ways. In Ref. |11|, use is made of the techniques described in 
Jl2| to obtain the kinetic equations for non-relativistic Wigner distribution functions of neutrinos. Alternatively, in 
references jl3j , |T4| it is derived the time evolution of a neutrino density matrix p constructed as macroscopic averages 
of generalized occupation numbers, in the way Pij(p) =< Pij(p) > , where Pij(p) = a^{p)ai{p) and aj~(p) (a,i(p)) is 
the creation (destruction) operator of neutrinos with flavor j (i) and momentum p. Here, one starts directly from 
Heisenberg's equation 

id t p=[p,H] (1) 



1 In this work we will concentrate on medium effects in neutrino propagation, and will not discuss production and absortion 
mechanisms. For these processes, the reader may want to consult 3|. 
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with a Hamiltonian H = Hq + Hint-, where Hq is the free Hamiltonian and Hi nt is the interaction piece. Equation 
(|l|) is then expanded perturbatively (after macroscopic averaging). With this at hand, the authors have studied the 
possibility of flavor conversion in a supernova core (see the references above for more details). 

Both procedures give rise to an expansion in powers of the Fermi's coupling constant Gf- The first term in this 
expansion (proportional to Gf) contains the modifications to the mass matrix due to the interaction, while the second, 
G F term, is the generalization of the Boltzmann collision integral to the case of flavor mixing. However, because both 
methods are based on non-covariant techniques, one does not obtain rclativistic dispersion relations for the neutrinos. 

In this paper, we make use of relativistic Wigner functions |l5f| , to describe propagation of neutrinos with flavor 
mixing in dense media. Relativistic Wigner functions have been successfully used to account for finite density and 
temperature effects in nuclear matter |l7J , jl8| , |nj . They provide an alternative to Green functions in a way which 
is well adopted to the development of kinetic equations. In addition to this, temperature effects are incorporated in 
a unique way, contrarily to the situation of Green functions. We also need to introduce correlation functions among 
neutrinos and the background (which, for simplicity, is considered to consist only on electrons). These correlations, as 
will be discussed later, take into account for the residual interaction of neutrinos and electrons beyond the mean-field 
approximation. We will obtain, in the next section, the equations of motion for these functions, under the assumption 
that the interacting background is in an equilibrium state. 

Next, we will consider some particular situations. In section 3, we examine the case when neutrinos are in equilibrium 
and correlations are neglected. The dispersion relations reproduce, in this case, the expected effective masses and 
in-matter mixing angles. This result supports our statement that neglecting correlations is equivalent to a mean-field 
treatment of the surrounding matter. 

In section 4 we examine some departure from the above simplest situation, by keeping spatial and/or time variations 
in the kinetic equations (with correlations still neglected) , and we consider propagation on a density- varying medium. 
In section 5 we examine with more detail a particular case, corresponding to a small effective potential and macroscopic 
inhomogeneities which are large enough. This is the situation encountered when dealing with solar neutrinos. By 
making the appropriate approximations, we recover the well-known formulae for the MSW effect. This shows the 
ability of Wigner function techniques to correctly reproduce both relativistic dispersion equations and transport 
equations on the same foot, at least in the cases we have studied. We end in section 6 by summarizing our main 
results and making some remarks. Some auxiliary results will be given in the appendixes. The construction of the 
Wigner function in the case of free neutrinos is showed in Appendix A. In Appendix B we analyze the neutrino 
dispersion relations and effective masses which appear within the Hartree approximation, and we extend the results 
of Appendix A to the construction of the corresponding Wigner function. 

In this work the metric is g^ v — diag(l, — I, — I, — I). We adopt the chiral representation for gamma matrices, and 
natural units {% = c = 1) are used. 

II. EQUATIONS FOR WIGNER FUNCTIONS 

In this section we derive the equations of motion for neutrino Wigner functions. Neutrinos are assumed to propagate 
on a matter background in equilibrium, an hypothesis which holds in all the astrophysical scenarios mentioned in the 
introduction. 

In order so simplify the equations as much as possible, we consider only two neutrino flavors (namely electron 
and muon neutrinos) . Our treatment can be generalized in a straightforward way to include more neutrino flavors. 
Also, our main concern in neutrino propagation is neutrino oscillations and flavor conversion, which are given by 
charged-current interactions with electrons. Therefore, we consider a background consisting on electrons, although we 
take into account both neutral and charged currents. As before, the formalism can be trivially extended to account 
for a more general situation, such as neutral-current interactions on protons and neutrons. In this paper, we treat 
neutrinos as massive Dirac particles in the most simple model for massive neutrinos, i.e. we treat them the same way 
as the all other fermions (leptons and quarks). Within this minimal extension of The Standard Electroweak Theory, 
the charge of the neutrino field which is conserved is the total lepton number L = L e + L^. On the other hand, as we 
deal with low-energy neutrinos (with energies of the order a few MeV), we adopt an effective contact interaction of 
neutrinos with the matter background. 

In what follows, neutrino magnitudes without a prime will indicate flavor states, and primes will be used for free 
mass eigenstates. In the next sections and in the appendixes, a tilde will be used for interacting eigenstates. Flavors, 
as well as mass eigenstates, will be labeled by lattin superindexes such as a and b. Spin subindexes will be generally 
omitted ; when needed, we use indices such as i, j, fc, ... to label them. Lorentz indices will be labeled by Greek letters. 
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Since we deal with two neutrino species, it is convenient to introduce vectors and matrices in flavor (or mass) space. 
We therefore define the neutrino and antineutrino vector fields^ : 

v & = ( V$) ) ' d( - x) s ( r & r & ) (2) 

We also introduce the following matrices in flavor space 

a-*- (3) 
n~ -(f|J) (4, 

With the notations A' 1 = 7 M (gy — <M7 5 ) and w' 1 = ^^(gv ~ 9A1 5 )- Here, the constant = — | + 2 sin 2 6*^ (5,4 = — | 

) corresponds to the vector (axial) contribution of weak neutral currents, 9w is the Weinberg's angle, while the 
constants gy = \ + 2 sin 2 6yy and 'gA = \ arise from neutral plus charged currents. 
With these notations, the Lagrangian density is written as : 

C(x) — C e (x) + C v {x) + Cj{x) 
with C e (x) the Lagrangian of free electrons, 

C v {x) = V(x)ij ll d ti i'(x) — V(x)Mv(x) (5) 
the corresponding Lagrangian of free neutrinos, and 

^ G G 

Ci(x) = fv(x)fl IJ -u(x)e(x)uJ tJ .e(x) pV(x)A^u(x)e(x)X lJ ,e(x) (6) 

V 2 v 2 

the interaction piece. In Eq. (||) M is the free mass matrix in flavor space [J From the above Lagrangian one 
readily obtains the equations of motion for the neutrinos : 

i'fduMix) - Mv{x) - ^n^d{x)%{x)u JfJ e{x) - ^A^(x)%{x)X ll e{x) = (7) 
v2 v2 



S„F(aj)*y + 9(a;)M + ^V{x)Q^{x)u^e{x) + ^9(x)A^t(a;)A p e(a;) = (8) 
v2 v2 



We now introduce the neutrino Wigner operator 

F^ v \x,p) = (2vr)- 4 | d*ye~™v){x + y/2)d?{x-y/2) (9) 
and the electron Wigner operator 

Fjf>(x,p) = (2tt)- 4 f cfiye-^ix + yfiySiix-y/i) (10) 
FromEq. (§) one readily gets that the Hermitian conjugate is given by 

F^\x,p)= 1 %Fl< v \x,p) 1 % (11) 



2 The symbol " on top of a magnitude means that we are dealing with a quantum operator. This will be used to distinguish 
this magnitudes from statisical averages. 
3 If flavor mixing exists, M is non-diagonal, its eigenvalues being the masses of free mass eigenstates. 
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Here and hereafter summation over repeated indices is understood. With the help of Eqs. (0) and (g) one can 
obtain the equations of motion obeyed by the neutrino Wigner operators. After some algebra, we get : 
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dF^(x,p) - 2ipF {lJ \x,p)\ + 2iMF ( - v \x,p) 



(12) 



-2iG I 

~7T 



(2tt)- 4 / d 4 yd 4 k [fle(y)we(i/) + Ae(y)\e(y)] F^(x,p- k/2)e lk ^ x ~> 



dF [u) {x,p) + 2ipF {,j) {x,p) \ 7 - 2iMF^(x,p) = 

(2tt)- 4 f d 4 yd 4 kF^(x,p- k/2) [QE{y)ue(y) + A?(t/)Ae(t/)] e - lk{y - x) 



(13) 
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"75" 



One could also derive the corresponding equations for the electron Wigner operator F^ (x,p). However, we will not 
need these equations under the approximations discussed in this paper, and therefore we will omit them. Of course, 
if one wants to investigate the next order to this approximation, the whole system of equations has to be taken into 
account. 

We are now interested in introducing statistical averages from the quantum operators introduced above. These 
statistical averages are called Wigner functions po| , and are the analogous to the distribution functions we need to 
describe many-particles systems. These are, in general, complex functions, and also contain a Lorentz structure which 
will be discussed later. The electron and neutrino Wigner functions are defined, respectively, as : 



F^(x,p) =<F^(x >P ) >= (2tt)- 4 



<iV"' w < %{x + y/2)e t (x - y/2) > 



(14) 



Ff"\x,p) ^< F;p"'(x,p) >= (2^ 4 y d*ye-™ <W 3 (x + y/2)l>:(x~y/2) > (15) 

Here, the symbol < A > means the average of a given quantum operator A over a basis of quantum states which are 
compatible with the macroscopical knowledge of the system. The latter determines a given density matrix operator 
p . Thus the averaging is performed according to 

< A >= Sp{pA} (16) 

where Sp means the trace performed over the quantum basis. By taking this average on Eqs. (|l|) and @ one 
arrives to the following equations : 



[j(d - 2ip) + 2iM] FM(jc,p) = - f_"( 4 I d 4 yd 4 kd 4 k'e lk{y ~ x) 



F^\x,p) [ 1 {d + 2ip) - 2iM] = ^ I d 4 yd 4 kd 4 k!e- ' lkl -y-^ 



iG F V 2 

[< nTr(ujF {e \y,k'))F {u \x,p-k/2) > + < A7V(AF (e) (y,k'))F^ (x,p - k/2) >] (17) 

iG F y/2 
W) 

[< F ( - 1, \x,p~k/2)nTr(LuF (e \y,k')) > + < F^(x,p - k/2)ATr(XF^(y, k')) >] (18) 
In the latter equations, the symbol Tr means the trace in spin indices. 

Let us now introduce the electron-neutrino A^^and neutrino-electron B^ ev ' correlation functions : 

4$ a \x,x',p,p>) ee < F^ ab (x,p)F£\x',p<) > ~F^ a \x,p)F^(x',p') (19) 
Bl;: )ab (x,x>,p,p>) ee < F^(x,p)F^ ab (x',p') > -F^(x,p)F^ ab (x',p') (20) 

Then Eqs. (0) and ^ can be rewritten as 

[7(0 - 2ip) + 2iM] F {u) (x, p) = - lGF ^} I d 4 yd 4 kd 4 k'e lk{y - x) 

(27T) 4 J 

[ATr(\B{y,x,k',p- k/2)) +Q,Tr(cjB(y,x,k',p- k/2)) + 

ATr{\F ( - e \y,k'))F^(x,p-k/2) + riTr{ujF ( - e \y 7 k / ))F ( - u \x,p~k/2)} (21) 
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FM(x,p) h{d + 2ip) - 2iM] = [ d^kd^k'e'^-^ 

(27Tj 4 J 

[ATr(A(x, y,p — k/2, k')X) + flTr (A(x, y,p-k/2, k')u) + 
AF {v \x,p- k/2)Tr{\F {e \y, k')) + VLF {v \x,p- fc/2)2V(wF< e > {y, k'))] 



(22) 



The two-point correlation functions defined above are not independent. In fact, one can prove that they are related 
by: 



A \jki ab *( X ' X '>P'P') = lip l°r B pqrs a { x ' i x iP iP)lqkl°st 

which can be written, in a short way, as^| 

A^(x,x',p,p')=B^(x',x,p\p) 
This implies that Eq. (E2) is actually the Hermitian conjugate of Eq. (pl|). 



(23) 
(24) 



III. SYSTEM IN EQUILIBRIUM. HARTREE APPROXIMATION 

In order to obtain some insight into the physical meaning of the neutrino Wigner function we will, in this section, 
investigate the situation when both the electrons and the neutrinos are in equilibrium, which we characterize by all 
statistical magnitudes as being time-space translationally invariant. This means that one-point functions can not 
depend on x , while two-point correlation functions can only depend on the difference of coordinates, i. e. we assume 
that : 



and 



F {e \x,p) = F {e \p) 
F {lJ \x,p) = F {lJ \p) 

A^ (x, x' , p, p') = A^ (x -x',p, p') 
B^ (x, x' , p, p') = B^ {x -x',p,p') 

By defining the Fourier transform 

A^ ve \k,p,p') = (2ir)- 4 J d A xe- lkx A {ve \x,p,p') 
(analogously for B tyev ^\ we obtain the equilibrium equations for the neutrino Wigner function 



PI 

Gf 
V2 



M--j=J d 4 k' (ATr (\F^(k'fj + QTr (u^F^(k'] 



FM(p) 



d 4 kd A k' 



ATr [\B (ev) (k, k',p+ k/2)) + flTr [uB (w) (k, k',p + k/2)) 



(25) 
(26) 



(27) 



(28) 



(29) 



F [v \p) 
Gf 



P- M y§ / 



V2 



d^kd^k' 



k' for (xF^ik')) A + Tr (^( e) (fc')) fi)) 
ATr (A {ve ^ (k,p + k/2, k')X\ + QTr (A {ve \k,p + k/2, k')ui 



(30) 



As before, Eq. ( |30| ) turns out to be the Hermitian conjugate of Eq. (p9|). This set of equations is obviously not 
complete, and one should add the equations which are satisfied by the correlation functions in looking for such a 



4 For a matrix having both spin indices and generation indices, its Hermitian conjugate is obtained by interchanging the 
generation indices and then taking the Hermitian conjugate in spin space. 
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complete set. However, in doing so there automatically appear three-point correlation functions. This procedure can 
be infinitely continued, so that one obtains, instead of a closed set, an infinite hierarchy similar to the BBGKY (after 
Bogoliubov-Born- Green- Kirkwood-Yvon) hierarchy of classical systems ^J, |2^]. For classical systems, one usually 
truncates this infinite chain by neglecting correlations of order higher than a given one, usually by showing that 
higher orders correspond to more rapid variations in space and time. The next step consists then in incorporating 
perturabatively the next-order correlations. We will use here the analogy with the classical situation, and will first 
examine the situation at the lowest order, i. e. when all kind of correlations are neglected. Such approximation is 
commonly referred to as the Hartree approximation. As we will show, the neutrino dispersion relations which arise 
from this approximation correspond to modifications of the neutrino propagator at the one-loop level fllCfl . 

The Wigner function of electrons in equilibrium can be calculated using standard techniques. Following [j23f one 
has 

F^(p) = (2n)- 3 5(p 2 - ml) [e(p°)f+(p) + 6(-p°)f-(p)} ( 7 p + m e ) (31) 

with m e the electron mass, p° the time-like component of the electron four-momentum p and 9(x) the step function. 
The functions f^~(p), feip) are the Fcrmi-Dirac occupation numbers of electrons and positrons, respectively. In the 
frame where the matter fluid is at rest, they read as0 : 

ft(P)= eM E^ ) + 1 (32) 



(in Appendix A a similar calculation is shown for neutrinos) , where E = ^Jp 1 + m 2 , p e is the electron chemical 
potential and (3 the inverse temperature (we set the Boltzmann constant ks — 1). We can now calculate the traces 



and integrals appearing in Eqs. (29),(|30|). After some algebra, it is easily obtained 

Tr ( d 4 kX^F (e Uk) - { TrJd 4 k\°F^{k)=g v n , . 

Jr i W -\Trfd i kXF^(k) = 0; i = 1,2,3 ( ™> 



Analogously 



Let us now define the matrix (in flavor space) 



Tr J d A kuj Q F {e \k) =g v n 
Tr J d 4 k^F^ (k) = 0; i = 1, 2, 3 



Tr fh^Hk) = „_ r \ ' ', "^,,y> "/I" n , (34) 



here, V n = V2G F g v n is the effective potential for neutral currents, and V = V + V n , where V = \/2GFn is the 
corresponding potential for charged currents, with n = Aj d A k(2n)~^8(k 2 — ml)(9(p°)f l t(k) + 9(—p°)f~(k))k° = 
n e — n s the electron (minus positron) number density. With these notations, and neglecting correlations, Eq. j29| ) 
can be cast under the form : 

[ 7p _M- 7 °i(l- 7 5 )$]F(p) =0 (36) 

(since electrons have been integrated out, the neutrino superscript in Wigner functions will be omitted in what 
follows, in order to make notations simpler). It is easily recognized in Eq. ([56|) the appearance of the left-handed 
chirality projector P L = — 7 5 ), as a consequence of left-handed interactions. Let us also introduce the right-handed 
projector Pjj = |(1 + 7 5 ). With the help of these two projectors, one can define the following components of the 
neutrino Wigner function : 

F L = PlFPr 
F R , - P R FP L 
F rl = PrFPr 

F L r - PlFPl (37) 



5 Making the hypothesis that matter is at rest is equivalent to consider a particular Lorentz frame such that the fluid four- 
velocity is u M = (1, 0, 0, 0). Results in other frames can be obtained by restoring the four-velocity u^, as discussed in JLq. 



G 



If we apply Pl and Pr on the left and right of Eq. (|36|) we arrive to the set of equations 

TpFjn - MF L = 
1PF LR -MF r - 7 °$F Lfi - 
1P F R - MF LR = 
1P F L - MF RL - 7 °$F L = (38) 

By combining the above equations, one finally arrives to : 

(p 2 - M 2 - 1P1 °<t>)F L (p) = (39) 



(p 2 - M 2 - 1 ®p 1 M<S>M- 1 )F R (p) = (40) 

An important remark must be made. We are here considering a hypothetical situation where neutrinos had time 
enough to equilibrate with the background matter. Under this assumption, right-handed neutrinos can be produced by 
different mechanisms, such as spin-flip or pair production. However, production rates are suppressed by a factor rn v jE 
, where E is the neutrino energy and m„ its mass. In the astrophysical scenarios we are considering, the production 
rate of right-handed neutrinos is small, and therefore they can be neglected. By this reason, we will concentrate on 
the left-handed component Ft, (p) . Consistently with this approximation, neutrinos will be treated under the extreme 
relativistic limit m v jE << 1. This will imply that the neutrino field can be considered, approximately, as consisting 
on negative-helicity neutrinos and positive-hclicity antincutrinos. For Wigner functions, this is shown in Appendix 
B, where the Wigner function will be explicitly calculated in the interacting case. 

The dispersion relation obtained from Eq. ( ]39| ) can be diagonalized, and one obtains the well-known expressions 
for masses and mixing angles in matter : this is also studied in Appendix B. 



IV. NON-EQUILIBRIUM SYSTEM. TRANSPORT EQUATION. 

In this section, we investigate the evolution of neutrinos when deviations from equilibrium situations arise. More 
precisely, we will consider that neutrinos are created and propagate through the matter background. This implies 
that Eqs. (SM) - (^7|) will not be imposed, and therefore time and spatial variations have to be considered. This 
represents an additional difficulty in solving the system of equations for Wigner functions and correlation functions. 
For the moment, we will only consider a simple case, where correlations are neglected. This will serve us to investigate 
the possibilities of Wigner function techniques in deriving neutrino transport equations, and will allow in the future 
to study more complicated frameworks. As we will see in this section, the equations arising in this context are 
appropriate to deal with neutrino propagation and flavor conversion in the Sun. 

We return to Eq. (|l|), and assume that electrons can be locally characterized by their temperature and chemical 
potential, in such a way that Eq. ( |3l| ) is still valid. 

By performing the same procedure as in Eq. (|39|), we can derive an equation for Fl(x,p), which is now 

[□ - 4(p 2 - A/ 2 ) - Aip^ + 2i$(x)7 AI 7°a M + 4$0)7%7° + (d^(x)) -f°]F L {x.p) = (41) 

The next step is achieved by decomposing the complete Wigner function F(x.p) into the Dirac algebra. By using 
the chirality projectors, as in Eq. ( p7| ) one can write Fi,(x,p) under the form 

F L (x,p) = ~(l- l5 )f L t l (x,p) 1 » (42) 

where fh^{x,p) is a matrix in flavor space[| and transforms as a Lorentz vector. If the left-handed projector Pl is 
used again, we see that fL^{x,p) can be expressed as 

JTW) = \Tr[Fl\x,p)^] = ^(2^)- 4 | d^ye-™ <d b L (x + ±y)m(x-\y) > (43) 



We notice from here that Jlij.{x,p) is an hermitian matrix. 
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In the latter equation, vl is the left-handed component of the neutrino held. We next analyze the equation of 
motion obeyed by fL^(x,p). In order to do this, we substitute Eq. (|4|) into Eq. ([HI). After some algebra, we obtain 
the equations 

- -Uf L + (p 2 - M 2 )f L - (p a f L +p- f L ) 
+ip^f° L -~ ($(i)/£) + '-V -h) =0 (44) 

- -Uf L + (p 2 - M 2 )f L - *{x) (pHl +pfZ) + x ($(x)/ L ) 

+ ip^f L -~ (*(aO/L) + (*(»)/£) + i$(x)px /l = (45) 

Eqs. (Q) and (45) are the basic transport equations to be solved on a general situation, with the help of appropriate 
boundary conditions. We will investigate the consequences of this set of equations in a future work. For the moment, 
as a test, we will show that, under the circumstances usually considered when the MSW is studied, we reproduce the 
known equations for this effect. 



V. MSW EFFECT 

We consider neutrinos moving along a straight line (for example, the radial direction of the star). According to 
this, we assume that is parallel to p. This allows us to write 

f L (x,p)=pf(x,p) (46) 
where f(x,p) is a new function. We also introduce, for convenience, 

f L (x,p) = \p\g(x,p) (47) 
Next, we assume that neutrinos are ultrarelativistic, and that the effective potentials V and V n in Eq.(j3^) satisfy 

V, V n <C p° ~ \p\ ~ lMeV (48) 

The last condition will concern the characteristic scale of spatial and time variations of the neutrino distribution 
function. This scale has a macroscopic size R, at least of the order of 1 Km, or even more (in the case of the resonant 
zone in the Sun, for example). In this case, we can make the following estimate : 

-jf - |V/S| ~ ^ (49) 

When R = 1 Km, then 1/R ~ 10~ 16 MeV. If we are interested on variations of the distribution function on the 
scale R, then, together with hypothesis Eq. (|Is| ) we can simplify Eqs. (ffil ) and (45) to give 

[p 2 - M 2 - (p° + \p\Mx)] (/ + g) + ip^(f + g)=0 

[p 2 - M 2 - (p° - \p\Mx)] (/ - g) + ip^(f ~g)=Q (50) 

In the latter equations, f(x,p) and g(x,p) are Hermitian matrices (in flavor space). We can perform a local transfor- 
mation for each of them in such a way that both become diagonal (and therefore real). In this way, one can easily 
check that two possibilities are open for the system Eq. (^) to have non-trivial solutions. The two possibilities are 

det [p 2 - M 2 - (p° ± \p\)$(x)] = (51) 

We recognize in Eq. (|l]) the dispersion relations for neutrinos and antineutrinos, as described in Appendix B within 
the Hartree equilibrium hypothesis. However, quantities in the latter equation depend on the coordinate x . This 
means that the neutrino mass eigenstates can be obtained locally from the Hartree approximation (which gives the 
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same as the MSW effect). We then have, for neutrinos, det [p 2 — M 2 — (p° + \p\)<S>(x)] = 0. This implies the condition 
f{x,p) = g(x,p). Therefore, = \p\f(x,p) ~ p°f(x,p) , and we conclude that 

mx,p)~p»f(x,p) (52) 

within the same approximation. Finally, we have for Eq. ([42] ) 

F L (x,p) = - To)p^f(x,P) (53) 

where f(x,p) , for ultra-relativistic neutrinos, obeys the following equation of motion 

[p 2 - M 2 - 2<S>(x)p°] f + ip^f = (54) 

We have discussed above the possibility of making a local transformation which diagonalizes f(x,p). We have 
exploited the fact that, under these circumstances, it becomes a real matrix. However, from the physical point of 
view, it is more convenient to introduce a different local transformation, in such a way that the factor inside the 
brackets in the latter equation becomes diagonal, i.e. we consider an unitary transformation given by the matrix 

U M (x) 

7= U\JU M (55) 

such that 

Ul(x)[(p 2 — M 2 ) — 2$>{x)p°]Um{x) =p 2 - M 2 (x) (56) 
where M 2 (x) = u\ { (x) \M 2 + 2&(x)p°] Um{x) is a diagonal matrix, which contains the local mass eigenvalues []: 

^-(^ i4>) 

On the other hand, one finds 

Ul(x) (id^f(x,p)) U M (x) = id»pj(x,p) - [f(x,p), Ul(x)id^U M (x)} (58) 
Thus the equation of motion reads 

id»pj(x,p) + (p 2 ~ M 2 {x))f{x,p) - i[f(x,p), uI(x) P ^Um(x)} = (59) 
Let us write explicitly the matrices f(x,p) and Um{x) by defining 

u M (x)^( cosd /f] -si*M*n (60) 

y ' ysmduKx) cos9m(x) J y ' 

f(xp)=(^ {X ' P) J ^ P) ) (61) 
I[X > P) - \j 2 \x,p) f 22 (x,p)J [bL) 

where the functions / 11 (a:,p) and f 22 (x, p) are real, while f_(x, p) and f 21 (x,p) are the complex conjugate of one 



another. After substituting Eqs. (H),© and (|1J) into Eq. (§|), one obtains 



ip^/ u (z,p) + (p 2 - Mlix)) r{x,p) - ipfd^Mix) (f 12 (x,p) + f 21 {x,p)) = (62) 



7 We use a tilde to represent magnitudes in the interacting eigenstates basis, as mentioned in Section 2. 
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rp 



l dJ 12 (x,p) + (p 2 - Mi(x)) f 12 (x,p) + iffd^Mix) (f n (x,p) - f 22 (x,p)) = 



(63) 



ip*d»f 21 {x,p) + [f - M 2 {x)) f 21 (x,p) + ip^d^Mix) (f l \x,p) - f 22 (x,p)) = 



(64) 



ip^d^f 22 (x,p) + {p 2 - M 2 (x)) f 22 (x,p) + ip»d^ M {x) (f 12 (x,p) + f 21 (x,p)) = 



(65) 



A remark is in order. As can be seen from the above equations, it is not possible, in a general situation out of 
equilibrium, to perform a transformation that makes both the neutrino Wigner function and the mass matrix diagonal. 
However, local mass eigenstates can be used as a useful physical basis to simplify the equations, as done in this section. 

We can write the latter system of equations in a more familiar way. By taking the real part on the first and last 
one, we readily arrive to the conditions 



p 2 -M 2 (x)\ f n (x,p)=0 
p 2 -M 2 (x)} f 22 (x,p)=0 



(66) 



which imply that f (x,p) (f 22 (x,p)) is non-vanishing only when p 2 = Mf(x) ( p 2 = M|(a;)). Similarly, from the 
two remaining equations it can be easily deduced that the functions f 12 (x,p) and f 21 (x,p) have to vanish unless 



the condition p 2 = \ 
substitute 



M 2 (x)+M 2 (x) 



is fulfilled. Therefore, in the equations of motion for these functions we can 



p 2 -M 2 (x) 
P 2 -M 2 (x) 



-A(x) 
-U(x) 



(67) 



where A(x) = Mf (x) — Mf(x) is the in-medium neutrino mass difference. Furthermore, in the ultra-relativistic 
case, we can approximate the operator p^d^ in all the Eqs. (|62|)-(|65|) by 







p"d^\p\- + \p\- = \p\D 



dx 



(68) 



We have defined D — j% + -§^- After some manipulations, Eqs. (|62|)-(|6q) can be cast under the simple matrix form 



iDf(x,t) = [f(x,t),H(x)] 



Here, 



H(x) 



-i6' M (x) 



A(x) 

4p| '"Ml 



(69) 



(70) 



is the effective Hamiltonian in the mass eigenstates basis, and O'm(x) = -^9m(x). Eq. ( |69] ) can also be written in 
the flavor basis, if we undo the transformation introduced in Eq. (pq), with the definitions 



f(x,t) ^ ( f f Zf x ]t) ) = U M (x)f(x,t)Ul(x) 



After a straightforward calculation, we arrive to the equation 

iDf(x,t) = [f(x,t),H(x) 

with 



H(x) 



a{x) (3 
/3 — a(x) 



(71) 



(72) 



(73) 
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is the Hamiltonian in the flavor basis, and we introduced the notations 



a(x) = (A cos 26»- A{x))/A\p\ 

(3 = A sin26»/4|pl (74) 

where A(x) = 2|p|\/2GV n(x) is the induced mass due to charged currents, Ao = m| — m 2 is the vacuum neutrino 
mass difference and 9 the vacuum mixing angle. 

From the system of equations (|72]) , we first obtain that 

iD[f ee (x,t)+f^(x,t)} =0 (75) 

which implies that the total number of electron plus muon neutrinos is conserved during the propagation, as 
expected. We can then write f ee (x, t) + f^(x, t) — K , where K is constant during the propagation : DK = 0. By 
manipulating the equations, one can also derive a third-order differential equation for f ee (x 7 t), which reads as 

aD 3 f ee (x, t) - a'D 2 f ee {x, t) + 4a(a 2 + f3 2 )Df ee (x, t) - 4a' f3 2 (^f ee (x, t) - ^\ = (76) 

where a' = -j^ct. The latter equation coincides with the one derived by Mikheyev and Smirnov R] to describe the 
evolution of the survival probability of an electron neutrino in a non-constant density medium. This agrees with 
our interpretation of f ee (x,t) as giving the distribution function (proportional to the number density for a given 
momentum) of electron neutrinos, and analogously for / MM (x,i) as the distribution function of muon neutrinos. All 
known results for the MSW effect within the situation considered here can be reproduced from the above equations. 

Eq. (|7^) defines the evolution of flavor distribution functions, and can be compared to the corresponding results 
derived from other treatments. As we mentioned in the introduction, there are (to our knowledge) two different 
approaches to neutrino propagation in dense media which make use of some kind of distribution functions. Both 
methods are based on perturbation techniques, assuming that the Hamiltonian can be separated into two terms : 
H = Hq + Hi n t, with Hi n t considered as a small perturbation. They both arrive to an equation with a term which 
is second order in Gf and corresponds to a non-trivial (i.e., non-diagonal in flavor) Boltzmann collision integral 
of neutrinos interacting with other particles in the background. This second-order term is absent in our approach, 
at least when correlations are neglected. Within this approximation, and neglecting small derivative termsn, the 



evolution equations of Ref. [ JL 1[ | coincides with Eq. (72). A similar comparison can be made using the results derived 
in @, g§. 

VI. CONCLUSIONS 

In this paper we have studied propagation of two neutrino species in dense media with flavor oscillations. We used 
an scheme which is based on the introduction of relativistic Wigner functions and correlation functions. Our aim 
is to develop relativistic kinetic equations, and to analyze the possibility that such approach will correctly describe 
both the relativistic dispersion relations of in-medium neutrinos, as well as the appropriate evolution equations for 
neutrino distribution functions. This allows us to treat in-medium neutrino masses and kinetic equations on an equal 
foot. We considered the medium as consisting on electrons, although other constituents can be incorporated in a 
straightforward way. 

By writing the equations of motion for Wigner operators and taking statistical averages, one arrives to an infinite 
chain of equations of the BBGKY hierarchy-type. This infinite chain has to be broken at some level in order to obtain 
a solution, and we have considered here the lowest-level of approximation, which consists in neglecting correlations 
(the so-called Hartree approximation). We have first examined the situation in equilibrium, which reproduces the 
well-known results for relativistic dispersion relations jnj. Next, we assumed that the electron background is in 
equilibrium, although its density needs not to be constant, and neutrinos propagate out of equilibrium. This gives 
rise to a set of kinetic equations which have to be solved, in a general situation, by taking the appropriate boundary 
conditions. 



8 As discussed in , these terms should only kept in the equations when strong inhomogeneities of size L are present, such 
that \tf\L ~ 1. See also Q and @. 
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In order to obtain some insight into the above equations, we have considered with some detail the usual MSW 
scenario, in which neutrinos are relativistic and the scale of space-time variations in distribution functions have 
macroscopic values (such as flavor conversion in the Sun or newly-born neutron stars). In this case, our system of 
equations reduce to the results obtained by other authors with the use of perturbation techniques |pT|| , fl3f| when only 
the first-order correction is considered. However, none of these methods has been showed to incorporate the correct 
neutrino dispersion relations. 

Now the question which arises is wether the inclusion of correlations into our scheme will lead to the same results as 
in the previous references, also for the second-order terms. Indeed, this seems to be the case, since correlations turn 
out to be proportional to the coupling constant Gf, and substitution into the BBGKY hierarchy will give corrections 
of the order G F , as one can see from Eqs. ( ^T| ) and (p2|). A more detailed study of the kinetic equations developed 
here will elucidate this question and, perhaps, give rise to new phenomena in the physics of neutrinos in dense media. 
This will be the subject of a future work. 

APPENDIX A: FREE SYSTEM WITH ONE AND TWO GENERATIONS 

In this appendix we give explicit formulae for the non-interacting neutrino Wigner functions. We first consider the 
case with only one generation, and we will generalize these results to the case of two neutrino generations with a 
non-diagonal mass matrix. Neutrinos are assumed to be in equilibrium. Of course, equilibrium can not be reached 
in the non-interacting case, but one can imagine a situation where a very weak interaction is added in order for the 
system to reach equilibrium. This hypothetical interaction can then be turned off without changing the equilibrium 
properties of the system. 

We start with one generation of free neutrinos in equilibrium. The Wigner function then verifies the equations : 

(7P- m)F M (p) = 

F {u) {p)(-ip - m) = (Al) 

where to is the neutrino mass and F^ v '(j>) the equilibrium Wigner function. As in previous cases, the above 
equations are the Hermitian conjugate of one another. We next multiply the first one by (773 + to), which gives 

(p 2 -m 2 )F ( ^(p)=Q (A2) 

and this implies that the Wigner function vanishes whenever the condition p 2 — m 2 = is not fulfilled. We introduce 
the Grand Canonical density matrix operator 

~ =z -i e - H H- ll t ) , Z = Tre -/3(5-Mi) (A3) 

where H is the Hamiltonian, L the lepton number operator and [i the neutrino chemical potential. 
We obtain, after standard quantization of the neutrino field in the helicity-states basis: 

Fl;\p) = (2^)- 4 J d 4 ye-™ < D 3 (x + y/2)v t (x - y/2) >= 

= (2tt)- 4 f d 4 ye- ipy (2t:)- 3 f d 3 k ^ [< N x (k) > u^(k)u^{k)e lky - < N x (k) > vf(k)v^(k) e -' lky ] 

•* ^ A=± 

= (2n)- 3 f d 3 k [< N x (k) > u$(k)u$(k)5(p - k)- < fix(k) > vf(k)v^(k)6{p + k)] (A4) 

J A=± 

where N\(k) ( N\(k)) is the number operator for neutrinos (antineutrinos) with momentum k and helicity Af]. 
The notations u x (k), v x (k) 7 ... for spinors and (k),a\(k), ... for creation and destruction operators have their usual 
meaning. The above averages can be calculated using standard techniques, giving : 

< N x (k) >= Tr(pa{(k)a x (k)) = -^-L^— ^ = f(k) (A5) 



9 Operators will be considered in normal order. 
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<ftx(k)>=Tr(pP x (k)b x (k)) = 



f(k) 



(A6) 



with Ek = V to 1 + k 2 . Let us define the matrices 



E£.(fc)=u*(fc)u*(fc) 

Z\k) = -v$(-k)v?(-k) (A7) 
More explicitly, if we introduce a coordinate system in such a way that k = k(sm9 cos </>, sin# sin0, cos 9) one has 

1 I m —Ek — k 



m 

-Ek + k m 



(A8) 



(A9) 



£+(fc) - 



2E k \E k +k m I PW 



(A10) 



m E k + 

2E k \E k - k m 



p(k) 



where 



, _ , cos 2 | | sin 9 e 
- I \saL9e** sin 2 | 



and p(fc) = fj,(—k). By substituting into Eq. (A.4) one obtains, after some algebra 

F<">(p) = (2 7 r)- 3 2 J B p( 5(p 2 - A/ 2 ) £ [0(p°)/(p)S A (p) + 0(-p°)/»S A (p)] 

By performing the sum over helicities, we finally arrive to the expression 

F (u) {p) = {lP + m)f w (p) 
where the scalar function fw(p) is given by 

fwip) = ^Tr [F^(p)] = (2n)- 3 6(p 2 - m 2 ) [0(p°)/(p) + 9(-p°)f(p)] 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



The generalization of the above formulae to more than one neutrino flavor is done by using mass eigenstates as an 
intermediate step. As defined in Section 2, we use a prime to represent such states. In the case of non-interacting 
neutrinos, they arise from diagonalization of the free mass matrix M. The equation of motion for the neutrino Wigner 
function is now, in the flavor basis 



( 7 p- M)F {v) {p) = 
The free mass eigenstates are related to the flavor states through a matrix 



U = 



cos a — sin f 
sin 9 cos 9 



(A16) 



(A17) 



in the way 



V(x) = Uv'{x) 



(A18) 
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and one has 

M = UM'U^ 

where M' is the eigenvalues matrix 

M = diag(mi, ma) 

The Wigner function in the mass eigenstates basis, F'( u > (p) is related to F^ u \p) trough 

F {v \p) = UF' M {p)U ] 

and will obey the diagonal equation 

(p 2 -M l2 )F'^(p) = 
By quantizing the massive fields v la (x) (a = 1, 2 for two generations) we obtain the result 



(A19) 
(A20) 

(A21) 
(A22) 



F\ 



(p) = (2vr)- 4 / d 4 ye~™ < V°(x + y/2)u?(x - y/2) 



>- 



(2tt)~ 3 / d 3 k [< > uy x {k)u°;- x (k)d{p- k) 



< N x (k) > v]' x {kX' A (k)5(p + k)}S a , b a, b = 1, 2 



a, A 



(A23) 



with similar notations as before. The superscript a has been added to label different mass eigenstates. In this basis, 
then, the Wigner function is diagonal, as expected from Eq. ( A22 ). The Wigner functions appearing in the latter 
equation are given by 



F'W u (p) = ( 27 r)- 3 % 2 -m 2 ) [0(p°)/i(p)+0(-P O )/i(p)j (TP + mx) 
F'M 22 (p) = (2n)- 3 S(p 2 -m 2 2 ) [9{p°)f 2 (p) + 9(-p°)f 2 (p)] ( 7 p + m 2 ) 

and we introduced the notations 

Up) ' 



/»(P) = 



1 



1 



(A24) 

(A25) 
(A26) 



(a = 1 , 2) w ith Ep = \/m? a + p 2 . We notice that the chemical potential is the same for both generations of neutrinos 
in Eq. ( A25 ). This is so because, in our model, there is only a conserve d cur rent : the electron plus muon lepton 
number. A similar comment ca n be made regarding antineutrinos, Eq. (A26). The Wigner function in the flavor 
basis is then obtained from Eq. ( A21). 



APPENDIX B: NEUTRINO MASSES AND WIGNER FUNCTIONS IN THE HARTREE 

APPROXIMATION 



We now discuss with some detail the dispersion relations and neutrino masses which arise from the equations of 
motion of the Wigner function within the Hartree approximation. In order to simplify the notations as much as 
possible, and to obtain insight into the problem as well, we start with one neutrino flavor (let us consider electron 
neutrinos). Then, Eq. (36) takes the simpler form 



[7p-r7i-7°i(l-7 5 )T/] J F l(l/) (p) =0 



(Bl) 



Dispersion relations are the necessary conditions for Eq. (Bl) to have solutions other than the trivial one F^(p) = 
0. One can compute the determinant for this equation by setting p = (po, 0, 0,p z )- The result then is 



det[7p-m-7°i(l-7 5 )V r ] = 



Po - v\ ~ m 2 - (p z + po) V 



Pi 



p 2 z - m 2 + (p z - po) V 



(B2) 
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By equating to zero this determinant, and solving for pq, one obtains four solutions : po = E + , po = £L, po 
and pq — —E-, where 



V 
1 



+ V(? ~Pz) 2 + m 2 



-E 4 



(B3) 



V / V 
E- = - + W(-+p z ) 2 + m 2 



(B4) 



£4 



V ,V 



Pz) 



(B5) 



E- 



V ,V 
-2 + V ( 2 



(B6) 



In order to identify the above results, we start from Eq. (Bl) and use the projection operators to obtain an equation 
for F^ l/ \p), as in section 3. One then arrives to 



(p 2 - m 2 - jjrf°V)Fl v> (p) = 
We take into account the following identity 



1P1 



p° +p - a 

p° — p ■ a 



This means that the positive-helicity component F~£(p) has to obey the relationship 

[p 2 -m 2 -V(p Q -\p\)}FZ(p) = 



(B7) 



(B8) 



(B9) 



The solutions for the dispersion relation p 2 — m 2 — V(po — \p\) = are given by po = E + and po = —E + . One can 
also make an analogous study for F£(p), which gives pa = E_ and po = —E^ as possible solutions. This allows us 
to interpret Eqs. ( |B3| )-(B6) as corresponding to neutrinos (E) and antineutrinos (E) of positive and negative helicity 
(subindexes + and — , respectively). One can check that these solutions, obtained from the Hartree approximation, 
coincide with the ones obtained from the poles of the neutrino propagator at the one-loop level (see, for example, 
@)- 

We now concentrate on ultra-relativistic neutrinos, and consider densities as occur in normal or compact stars. 
Under these conditions, we can neglect terms such as V| "p* | and m 2 by comparison to |~p*| 2 and, in this way, the 
dispersion relations give us the neutrinos effective masses for each degree of freedom, which are 



Mi 



(BIO) 



M 2 _ =m 2 + 2V\^ 



M + =m 2 - 2V\-f 



(Bll) 
(B12) 



M 2 _ = m 2 



(B13) 



We used the same notations as for the energies. In these equations we see that positive polarization neutrinos and 
negative polarization antineutrinos behave, approximately, as free particles. 



15 



The above results can be generalized to the case of two neutrino flavors in a straightforward way. Let us start from 
Eq. (^) and transform it to the mass eigenstates basis in vacuum, as discussed in Appendix A for the non-interacting 
case. We then arrive to the following equation 

{p 2 - M' 2 - 7 p7 $')^ (l/) (P) = (B14) 

Following the discussion we made above for one generation, we readily obtain the equation for the negative-helicity 
component (we replace po to £L for neutrinos) 

(p 2 - M' 2 - (£L + m)&)F£'M(p) = (B15) 
The latter equation suggests us to define an effective mass matrix M. _ as follows 

M 2 _ =M' 2 + (E_ + \0\)& (B16) 
If flavor mixing exists, this matrix is non-diagonal. More explicitly, we find 

U 2 _ { (V n + V cos 2 9)(E_ + \p\) + m 2 -V sin6 cos6(E_ + \p\) \ 

Jvl --y -V sin 6 cos 6(E. + \p\) (V n + V sin 2 0)(£L + \p\) + m 2 J {aL<) 
The corresponding dispersion relation is then 

det(i> 2 -.M 2 ) =0 (B18) 



Eqs. (B16) and (B18)can be used to obtain the exact energy levels of neutrinos (and antineutrinos) as a function 
of the neutrino momentum. In the ultra-relativistic limit, one can approximate 

A >2 ( A n V / Acos 2 6 + mj -Asmdcosd \ , . 

M -~{0 A n ) + {-Asin9co S Asm 2 9 + m 2 ) (B1J) 

here, A = 2\p\V is the induced squared mass due to charged currents and A „ = 2\p\V n the analogous magnitude for 
neutral currents. We first consider the first term on the second hand of Eq. ( BlSp in order to find the eigenvalues of 
M?_ (following the notations introduced in Section 2, we represent the diagonal form of A4- by M ). Since neutral 
currents are diagonal in flavor, they can be added at the end. After some trivial algebra, the mass eigenvalues M\ 
and Mi are given by 

M l,2 = \ ( A + s ) T 2 K A o cos 20 - A) 2 + Aq sin 2 26] 5 + A n (B20) 

where E = m 2 , + m 2 and Ao = m| — m 2 . One can also obtain, by diagonalizing M 2 _ the in-medium mixing angle 
6m- The well-known result 

sin 2 26 M = A ° Sln2 ^ - . (B21) 

(Ao cos 26 - A) 2 + A^ sin 29 V 1 

is then reproduced. 

The procedure discussed in Appendix A can be extended in order to construct the neutrino Wigner functions in 
the equilibrium state described by the Hartree approximation. We will only give the final result. As before, it is 
illustrative to consider first the case of one generation (for example, electron neutrinos). One obtains 

FM (p) = F- (p) +F+(p) + F- (p) + F+ (p) (B22) 

where 

F~(p) = (2n)- 3 2E p S(p 2 - M 2 _) f (p)S- (p) 
F+(p) = (2n)- 3 2E p S(p 2 - M 2 )/(p)E+(p) 
F~(p) = (2^2E p 8{p 2 - Ml)f(p)Z-<j?) 

F+(p) = (27Ty 3 2E p S(p 2 ~ M 2 + )f(p)t+{p) (B23) 
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and the effective masses are defined in Eqs. (B1C)-(B13). The chiral left-handed component F L u \p) , as defined in 
Eq. ( |37j ) can be approximated by : 

F«(p) ~ F£{p) + F+(p) = (2n)- 3 2E p [S(p 2 - M 2 _)f(j>)^(p) + 6(p 2 - M 2 )/>)£+(£)] (B24) 
We have introduced the following definitions : 

'0 



V- L {p) = P L Y,-{p)P R 
t+{p) = P L t+(p)P R 



1 



1 



pip 
pip) 



(B25) 



£+($ = P L X+(p)P R ~ 

The latter approximations hold for ultra-relativistic neutrinos. Thus, under this approximation the field contains only 
two degrees of freedom : neutrinos with negative helicity and antineutrinos with positive helicity. 

We now go to the case of two neutrino generations. As in the free case, the Wigner function will be diagonal if 
considered in the interaction eigenstates basis : 



F^ ab {p) 



_ ( F^ip) 











FH 22 (p) 



(B26) 



Moreover, each one of the diagonal components can be easily obtained by taking into account the corresponding 
dispersion relations. One then arrives to a set of equations similar to Eqs. (B22)-(B23) for each one of the two 
diagonal components of the Wigner function. As in the case of one generation, we concentrate on the chiral left- 
handed component, which is finally approximated by 



f l ip) 



' — ab 

Fl ip) 



Ft. 



\-ab 



(p) = {2iry 3 2E p 6 ab 



6( P 2 - M 2 a )f a {p)^- L (p) + 8{p 2 - M a )f a (p)E+(p) 



(B27) 



with a, b = 1,2 . The functions f a (p) and f a (p) are the same functions defined in Eqs. ( A25| ) and (A26), but with 
the replacement m a — > M a . The two terms in the latter equation obey the relations 



( 1P ~M ab )F^ ab (p) 

--^- ab ab 

( 1P -M )F L (p) 



(B28) 



where M ah and M are the (diagonal) effective masses matrices for neutrinos and antineutrinos. In order to construct 
the Wigner Function in the flavor space, we perform an unitary transformation Um that is defined by the mixing 
angle in mater 9m ■ In this way we have, for example : 



\V) 



cos 2 (9 M )Fr 11 



ip) 



s\n 2 (6 M )F7 22 



F-^(p) = F-» e (p) = sm(0 M )cos(0M) 



F£ U 



(P) 
ip) 



ip) 



F~^(p) 



-11 



■>M 



)Fl (p) + ■ 



—22 

hi)F L (p) 



(B29) 
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